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■ We present a consistent theoretical approach for the study of nonequilibrium effects in chiral fluid 

dynamics within the framework of the linear sigma model with constituent quarks. Treating the 
quarks as an equilibrated heat bath we use the influence functional formalism to obtain a Langevin 
fT^ ■ equation for the sigma field. This allows us to calculate the explicit form of the damping coefficient 

and the noise correlators. For a selfconsistent derivation of both the dynamics of the sigma field and 
the quark fluid we have to employ the 2PI (two-particle irreducible) effective action formalism. The 
energy dissipation from the field to the fluid is treated in the exact formalism of the 2PI effective action 
where a conserved energy-momentum tensor can be constructed. We derive its form and comment on 
^ ' approximations generating additional terms in the energy-momentum balance of the entire system. 

I. INTRODUCTION 

> . 

, Since neutron stars are too far in space, the Big Bang is too long back in time and astrophysical observations are 

psj ' too indirect, the only possibility to investigate the phase transition of QCD experimentally is by heavy-ion collisions. 

\^ Here, matter is created under extreme conditions, in small systems and with fast dynamics. It is, therefore, likely 
that nonequilibrium effects play an important role in the evolution of the fireball. Nonequilibrium methods have 
also been applied to other systems of extreme conditions, e.g. in the inflationary evolution of the early universe 
and the formation of Bose-Einstein condensates |3, 4]. 

^""^ I The predictions of most observables of the conjectured critical point are based on the striking feature of diverging 
. fluctuations and correlation lengths in thermodynamic systems at criticality. These should be visible in event-by- 
^ event fluctuations of particle multiplicities [5, 6]. The proposed experimental measures of fluctuations, which are 
based on the second moments of the particle distributions did not show any increase in the energy range of E]^ = 

rS 18 — 158 GeV at NA49 [7]. The key observable for the beam energy scan at RHIC is the net-proton kurtosis fa,^, 
[ based on the fourth moment of the distribution. It has the advantage that in thermodynamic systems the kurtosis 
diverges with higher powers of the correlation length than the second moments |9]. It is, however, unclear if the 
system created in a heavy-ion collision is in thermal equilibrium at the phase transition. Here, relaxation times 
become large, whereas the dynamics of a heavy-ion collision is very fast. Even if the system is in equilibrium above 
the phase transition it is likely to be driven out of equilibrium as it cools through the phase transition. At a second 
order phase transition this is called critical slowing down and severely limits the growth of the correlation length [10] . 
As a consequence, one expects any signal of a critical point to be weakened in a dynamic nonequilibrium situation. 
Signals of a first order phase transition, however, are based on nonequilibrium effects, such as supercooling. Here, 
parts of the system remain for some time in the high-temperature phase even below the transition temperature. 
The decay of this unstable state occurs via either nucleation or spinodal decomposition. A proposed signal for a 
nonequilibrium situation at a chiral phase transition is the enhancement of soft pions from the decay of disoriented 
chiral condensates flTI - ITsll . 

Nonequilibrium phenomena can be studied in the framework of Langevin dynamics. Here, the chiral fields are 
split into hard and soft modes. The hard modes constitute a heat bath for the soft modes. Due to the interaction with 
the hard modes the soft modes undergo additional dissipative processes. Calculations of the influence functional 
for a system of soft modes interacting with a heat bath of hard modes have extensively been studied in cp"^ theory 
llT6l - [l9.] , in gauge theories l,20i,,2l,1 and in 0{N) chiral models [22]. However, the back reaction of the soft modes on 
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the hard modes, e.g. the exchange of energy is typically not included in these influence functional studies. 

A framework which is based on the conservation of energy and particle numbers is fluid djmamics. Chiral fluid 
dynamics combines the dynamics of heavy-ion collisions with the explicit propagation of the chiral fields at the chiral 
phase transition. The common idea is to embed a low energy effective model of QCD into a fluid dynamic simulation 
of a heavy-ion collision. A model that is particularly well suited for this combined approach is the linear sigma model 
with constituent quarks ll23ll . The dynamics of the quarks is reduced to a fluid dynamic evolution of densities. This 
gives rise to the name chiral fluid dynamics. The quarks are thus assumed to be locally equilibrated and act as a 
heat bath with local temperatures. In [24] an additional dilaton field was included in the explicit dynamics to model 
scale invariance breaking in QCD. Strong nonlinear oscillations were found for all the classical fields at the phase 
transition. From these oscillations the particle production of sigmas, pions and glueballs was calculated. In [2^ 
initial fluctuations were propagated deterministically through the phase transition. The linear sigma model with 
constituent quarks treats the quarks and antiquarks and the mesons on equal footing. In the real world confining 
forces recombrne quarks and antiquarks in mesons and baryons below the confinement critical temperature. The 
aspect of confinement is not included in the linear sigma model with constituent quarks. In extensions of the model 
gluons are included on the level of the dilaton field [24] or the Polyakov loop [26, 27]. In all the existing approaches 
so far ||2^|2^[28.] the time evolution of the chiral fields is given by the classical equations of motion. Thus, the chiral 
fields are explicitly propagated out of equilibrium with the quarks and do not relax to their equilibrium value for 
the temperatures given by the heat bath. Relaxational and stochastic processes have been neglected. 

A formalism which is capable of including dissipation and noise together with a consistent back reaction on the 
heat bath is provided by the so-called 2PI formalism. Two-particle irreducible (2PI) diagrams are closed diagrams 
which do not fall apart when cutting two lines. Consequently the 2PI formalism deals with full one-point (mean- 
field) and full two-point functions (full propagators). It has been originally developed as a functional formalism in 
II29II and was extended to yield transport equations in a quantum theory that preserve the macroscopic conservation 
laws and are thermodynamically consistent in the equilibrium limit [30-32]. These can be found from a certain class 
of approximations, called O-derivable, to the Schwinger-Dyson equation for the two-point function. In the path 
integral formulation of this approach 1 33|| it turns out that the O-f unctional is identical to the nontrivial part of the 
two-particle irreducible effective action. In [l33] this selfconsistent approach is generalized to arbitrary nonequilib- 
rium many-body systems. Renormalization of O-derivable approximations is studied in jssl - ls^ . 

In this paper we extend existing chiral fluid dynamic models by the consistent inclusion of damping and noise in 
the dynamics of the order parameter of chiral symmetry, the sigma field. To achieve this we apply two methods of 
nonequilibrium quantum field theory, the influence functional and the two-particle irreducible effective (2PI) action, 
to the linear sigma model with constituent quarks. In the influence functional method the separation of the systems 
seems obvious in our approach to chiral fluid dynamics. We consider the quarks as the heat bath and treat the 
sigma field as the relevant sector. The influence functional is obtained by perturbatively integrating out the quarks. 
We obtain the proper Langevin equation for the sigma field. However, in the influence functional method it is not 
possible to control the local equilibrium properties of the quarks without further assumptions. We need to go beyond 
existing studies of Langevin equations by putting special emphasis on the properties and the evolution of the heat 
bath itself. We expect that the back reaction to the heat bath induced by the dynamics of the chiral fields can be 
important for the overall evolution of the system. 

In order to derive the nonequilibrium propagation of the chiral fields and the thermodynamic properties of the 
quarks selfconsistently we apply the formalism of the 2PI effective action. The analysis of the chiral order parameter, 
the a field will be restricted to the mean field. One often defines the mean field as an average over quantum and 
thermal fluctuations in which case there is neither damping nor noise. Within the 2PI effective action formalism 
the mean field is obtained from an integration over quantum fluctuations only and thus still contains the necessary 
information about dissipation and noise. 

For a first qualitative analysis we restrict our model on the effect of the chiral phase transition and do not include 
the Polyakov loop. The Polyakov loop extended model describes confinement well on a statistical level and thermo- 
dynamic quantities have been calculated i26ll27ll but it is not evident how to treat the Polyakov loop dynamically. In 
our present analysis we also neglect the pion degrees of freedom. Their vacuum expectation value vanishes anjrway. 
From a conceptual point of view their inclusion provides no further problems. Their physical significance is largest 
below the chiral phase transition. As quasi-Goldstone bosons the pions are very light. Therefore they can be pro- 
duced from the decay of the sigma meson, leading to additional damping of the sigma field. The hard modes yield a 
sizable contribution to the heat bath with their own equation of state. Thus this needs to be included in a two-fluid 
dynamic description. Obtaining the Langevin equation for the soft pion modes is also a straightforward task. The 
inclusion of pions is devoted to future work. In the present work we explain the formalism of the coupled system 
of the explicitly propagated degrees of freedom and the fluid dynamic expansion of the heat bath. In an upcoming 
paper we will present detailed numerical results in a one-fluid dynamic set up. For the time being we concentrate 
on the heat bath of quarks and antiquarks. 
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This paper is organized as follows. In section|ll]we briefly review the linear sigma model with constituent quarks 
and its phase structure. The influence functional is calculated in section |III] including the damping and the noise 
kernel. In section |IV] we apply the formalism of the 2PI effective action. Both approaches yield the same Langevrn 
equation for the sigma field, calculated in section |yl In the exact formalism of the 2PI effective action a conserved 
energy-momentum tensor of the entire system can be constructed. We will derive this form and comment on the 
energy-momentum conservation of approximations to the full equations in section IVll 



II. THE LINEAR SIGMA MODEL WITH CONSTITUENT QUARKS 

The linear sigma model with nucleons has been studied for years as the prototype effective model of dynamic 
chiral symmetry breaking. The coupling to d5mamic quark degrees of freedom instead of nucleons yields the addi- 
tional feature that the quarks with light current masses obtain a heavy mass at the phase transition and thus turn 
into constituent quarks [38-41]. The Lagrangian reads 

C = -q[iY^f, - g{cr + ij5Tn)]q + ^{dficrdf'a) + ^{d^.Tidf tz) - U{cr, n) , (1) 

where q = {u, d) is the constituent quark field, a the sigma field and ft the pionic fields. The strength of the coupling 
between the quarks and the chiral fields is g. In the vertex for the pion-quark coupling the 75 matrix enters to account 
for the pseudoscalar nature of the n mesons and the isospin Pauli matrices t for the isospin degeneracy of the pions. 
The interaction between the chiral fields is given by the potential 

-,2 2 

U{a,Ti) = — [a^ + jf -v^') -h^a-Uo. (2) 

The Lagrangian ^ is invariant under SU(2)l x SU(2)r symmetry transformations if the explicit symmetry breaking 
term hq vanishes in the potential ||2]l. The parameters in ||2ll are chosen such that chiral symmetry is spontaneously 
broken in the vacuum, where {a) = = 93 MeV and {ft) =0. The explicit symmetry breaking term taking 
into account the finite quark masses is hq = fn^n^ with the pion mass niji = 138 MeV. With these requirements, 
1/^ = — in^ I A^. Choosing 7^ = 20 yields a realistic vacuum sigma mass = 2A^/^ + w 604 MeV. In order to 
have zero potential energy in the ground state the term LIq = m\l (4A^) — f^m^ is subtracted. At a coupling g = 3.3 
the constituent quark mass in vacuum is niq = 306.9 MeV. 

The one-loop thermodynamic potential in mean-field approximation is 

n(T,^) = -^lnZ= U((7,7r) + n,,-, (3) 
with the chiral potential U{(t, n) and the quark contribution 

n,,(T, ^) = -d, I ^ (^E + T In (1 + exp ) + T In (1 + exp (^^) ) ) , (4) 

where }i is the quark chemical potential and dq = 12 is the degeneracy factor of the quarks for Nf — 2 flavors, Nc = 3 
colors and the two spin states. Though technically very similar there is a noteworthy difference to the calculation of 
the partition function for free fermions with mass ntf. In the Lagrangian ([ij there is no fermionic mass. The quark 
mass is generated by nonvanishing expectation values of the chiral fields due to spontaneous symmetry breaking. 
During the evaluation of the functional determinant in Dirac and isospin space one generates a term defined as the 
effective mass of the quarks 

ml,, = g\cr^ + n^). (5) 
Then, the energy of the quarks and antiquarks is 

E = + ml^f = ^f + g^{(T^ + n^) . (6) 

This is obviously not a medium-independent quantity as the chiral expectation values depend on both medium pa- 
rameters T and Thus, the divergent term in l|4ll cannot be subtracted as a simple zero- temperature contribution. It 
needs to be renormalized more carefully. By using standard renormalization techniques one part of the divergence 
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can be absorbed in the parameters A and v of the classical potential U{cr,fz), while a logarithmic term depending 
on wZgff and the renormalization scale remains. In [42J it was shown that by neglecting this contribution one fails 
to reproduce the second order phase transition for p = in the chiral limit. In [43] the renormalization scale de- 
pendence was investigated phenomenologically. A thorough study of medium dependent corrections to mean-field 
calculations, perturbative and renormalization group approaches to Yukawa theory, also shows a crucial effect on 
the phase structure jiil-liill. 

To achieve the goal of this work, namely the coupling of chiral nonequilibrium d5mamics at the phase transition 
to a fluid dynamic expansion of the matter, we need a field-theoretical model exhibiting a phase transition. This is 
given by the mean-field approximation |47] and we can neglect the effects of the vacuum correction. 

The pressure of the system is 

p{T,ii) = -n{T,fi), (7) 

from which all thermodynamic quantities can be calculated. We are especially interested in the energy density. It is 
given by the thermodynamic relation 

e{T,ii) = Ts-p + }Ln, (8) 
with the entropy density s = {dp/dT)f, and the baryon density n = — (dp/dp) j. Then 

f^p{T,^) 



p{T,}i)-}i 



(9) 



The current understanding of the phase diagram of QCD is based on lattice results, which reliably predict that 
the phase transition is a crossover for zero baryochemical potential fi^, and on less settled model studies, which 
claim that at larger baryochemical potentials and lower temperatures the phase transition is of first order [49] . Then, 
this first order phase transition ends in a critical point. The phase structure of the linear sigma model with con- 
stituent quarks in mean-field approximation is qualitatively the same. For a first qualitative analysis, we can fix the 
baryochemical potential at ^( = and tune the strength of the phase transition by changing the coupling constant g 
iSOHSlfl. At /; = the thermodynamic potential is 

n{T,cp) = U{a,n)+n^^ = U{a,Tt)-2d^T J log (^1 + exp (^-|) ^ . (10) 

The sigma mass is given by the curvature of the thermodynamic potential at the equilibrium values of the chiral 
fields 

'"^ = |^l'-=-e,- (11) 

For the realistic coupling, g — 3.3, and ^ = the effective potential changes smoothly from the high-temperature 
phase to the low-temperature phase, see figure For higher couplings g the effective potential starts to exhibit a 
first order phase transition. In figure l|2]l we show the effective potential for various temperatures and g = 5.5. 

Above the critical temperature there is one minimum close to a 0, which is the minimum without explicit 
symmetry breaking. At a critical temperature of Tc = 123.27 MeV the two minima are degenerate and form the 
two coexisting phases. Below Tc the high-temperature minimum becomes unstable but exists down to the spinodal 
temperature T^p = 108 MeV jH. 

If one carefully looks for the vanishing of the barrier by decreasing the coupling g, the effective potential shows 
the shape of a second order phase transition for g = 3.63 with a critical temperature of Tc = 139.88 MeV. Here, the 
curvature at the minimum becomes very flat, see figure |3l and the sigma mass decreases to = 26.6 MeV. This 
value of nicr is significantly smaller than the vacuum mass, but still not zero. In order to lower the sigma mass further 
it would be necessary to time the parameters g and Tc more precisely. 

For all three phase transition scenarios we show the temperature dependence of the mass of the sigma field in 
equilibrium fWu in figure|4j 



III. THE INFLUENCE FUNCTIONAL FOR THE LINEAR SIGMA MODEL WITH CONSTITUENT QUARKS 

The influence functional method |52, 53] gives a reduced description of the entire system with focus on the evo- 
lution of the relevant variables. The details of the environment are eliminated by integrating out the environmental 



5 



T=160 MeV 




-0.4 -0.2 0.2 0.4 0.6 0.8 1 

a/T 



Figure 1. The effective potential for a coupling g = 3.3 and three different temperatures. The potential changes smoothly between 
the high and the low temperature phase. It describes a generic crossover transition. 
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Figure 2. The effective potential for a coupling g = 5.5 and temperatures T > Tc, T < Tc and the critical temperature. At the 
critical temperature the two minima are degenerate and represent the two coexisting phases for a first order phase transition. 



fields in a path integral over the closed time path contour fsi*, "SS*]. For the linear sigma model with constituent 
quarks we assume the following splitting: the irrelevant degrees of freedom are the quarks and antiquarks, which 
constitute the environment, or heat bath, and which are propagated fluid dynamically on the level of densities, and 
the relevant sector is that of the chiral fields, which we propagate explicitly. Again, we keep the pion fields fixed at 
their vacuum expectation value (n) =0. The free action of the quarks and antiquarks reads 

So[q,q] = i J d'^xq{x)'Yf'df^q{x) , (12) 

and the interaction between the quarks and the sigma field is of Yukawa type 



Smt[q,q,cr] = -g d^xq{x)q{x)a{x) 



(13) 
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Figure 3. The effective potential for a temperature above, a temperature below and the critical temperature at a coupling g = 3.63. 
At the critical temperature the minimum becomes very flat. This indicates a second order phase transition. 




Figure 4. The equilibrium value of the sigma mass nig- for the three scenarios: with a first order phase transition g = 5.5, a critical 
point g = 3.63 and a crossover g = 3.3. We also observe a discontinuity in 7Ti(j- at the first order phase transition. 



The influence fxinctional is calculated from a perturbative expansion in g. This might seem doubtful because g is 
of order 0{1). It must, however, be looked at the individual processes connected with the orders of the expansion. 
Higher orders in g involve more sigma modes and quark-antiquark pairs. The lower the density of the system the 
less likely these processes become and, thus, contribute less to the damping of the sigma field. Still, this is a crucial 
point in the application of the influence functional method. The explicit calculation of the influence functional to 
second order, which is the first nontrivial order, using standard techniques on the Keldysh-contour is deferred to 
appendix [a] It is convenient to write the influence functional in center a and relative variable Aa of the upper 
and the lower (— ) branch of the contour. It has one term that is linear and one term that is quadratic in Ac, 

iSu:[d;Acr] = i Jd'^xD{x)Aa{x) - i Jd'^x Jd*yA(T{x)J\f{x,y)Acr{y) , (14) 
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with the damping kernel 

D(x) = £Vy^(y) (^S< (x - y)S> [y - x) - S> {x - y)S< (y - x)^ , (15) 

and the noise kernel 

= |^s<(^ _ y)s>(j^ _ ^) + s>(x - y)S<(y - x)^ . (16) 

The influence functional 

Sif[^,A(7]= Jd^xDix)Aa{x) + ^JdL^xJd^yAo-ix)J^{x,y)Aa{y), (17) 

has an imaginary part. It is exactly this term that causes the underlying quantum system to decohere and allows for a 
classical description of the system's dynamics. This means that every trajectory can be assigned a unique probability. 
For these trajectories the relative field variable vanishes A(j(x) = 0, because trajectories that have a significantly large 
A(7 are exponentially suppressed. Together with the coarse graining of the environment, the decoherence leads to the 
quantum-to-classical transition of the system. Fluctuations in the classical equation of motion appear as a remnant 
of coarse graining and decoherence [56-60]. 

The zeroth order in the perturbative calculation of the influence functional gives one due to normalization and the 

(2) 

first order vanishes. Then, the second order contribution Sjp is equal to z'-times the influence functional itself, as can 
be seen by taking the logarithm of iA3i 

zSiF K, = ln(l + Sjp^' [a+, a-])- s\^^ [a+, cf-] . (18) 

We also note that in the perturative approach to the influence functional we neither obtain the thermal mass correc- 
tion nor the equilibrium properties of the heat bath. For the calculation of the mass correction one needs to include 
further information. For example, it is possible to find this term by directly calculating the equation of motion as it 
was done in |22|| for cp"^ theory. In the formalism of the 2PI effective action this is included automatically. 



A. The noise kernel and fluctuations 



At first glance, the noise kernel seems to be redundant, because it is quadratic in A(j(x) and, therefore, vanishes 
after variation with respect to A(r(x). However, the semiclassical concept of obtaining the equation of motion from 
varying the action is well-defined only for real actions. Here, the noise kernel introduces an imaginary part, which 
we need to rewrite in order to obtain a real action. This is done by introducing a new stochastic field ^, which 
discloses the physical significance of the noise kernel. The imaginary part of Sjp can be rewritten by making use of 
the Gauss integral evaluation 

expJ^-i^dS Jd'^yAa{x)Af{x,y)Aa{y)^ = j V^P[^]exp (^i Jd'^x^{x)Aa{x)^ . (19) 

Here, the stochastic weight P[^] is a Gauss distribution 

P[?] = N'exp (^-1 jd^x Jd\jaxW-Hx,y)ay)^ , (20) 

with a normalization constant N'. Then, the stochastic field ^ is fully determined by its first two moments, a vanish- 
ing expectation value and the variance: 

(?W> = 0, (21a) 
(ax)^(y))=Ar(x,y). (21b) 



This stochastic force f plays an essential role in the equilibration of the classical fields. By the dissipation-fluctuation 
theorem it enforces the relaxation to the correct equilibrium state II6III . 
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B. The semiclassical equations of motion 

The semiclassical equations of motion for the sigma field are obtained from the stochastic influence functional Sip, 
defined in 

exp(zSiF) = j V^P[^] exp jd^x{D{x) + ^(x))Ac7(x)^ = j V^P[^] exp (/Sip) , (22) 

by varying 

SclK]-Sdk1+SiF[^,Ac7] (23) 
with respect to AtJ and then setting Act = 0. From the classical action one obtains 



<5(SdK]-Sdk1) 



(24) 



5^a 

The semiclassical Langevin equation for the sigma field is 

-D = ^. (25) 



6a 

Note that the damping kernel D generally depends on a, too. 

IV. THE 2PI EFFECTIVE ACTION FOR THE LINEAR SIGMA MODEL WITH CONSTITUENT QUARKS 



The formalism of the 2PI effective action I29l - l32l1 is well-suited for our purpose because it yields a selfconsistent 
and thermodynamically consistent description of the entire system. Thus, it allows for a well defined back reaction 
of the relevant modes on the heat bath. In the scheme that is developed in the following we use the semiclassical 
approximation for the sigma field. We restrict ourselves to the sigma mean-field and do not include the propagator 
of the sigma field. We, thus, work with a theory of fermions coupled to an external mean field. Since the quarks have 
a vanishing mean field, they are represented by their propagators. Then the 2PI effective action is a functional of the 
sigma mean-field a'^{x) and the full quark propagator S'^^{x,y) 

F [a, S] = Sd [a] - fTr In S - iTiS^^S + F2 [a, S] , (26) 

where the trace operation includes Tr = d^x X^flavor I^Dirac ^^'^ ^cl W] is the classical action of the sigma mean-field. 
The free propagator for a fermion with mass m jr inverts the differential operator of the free theory 

{id-mf)S"J'{x,y) = -iSf{x-y), (27) 

by which it is defined up to the boundary conditions. 

The first three terms in | |26l l are the one-loop results. The additional term F2 [c, S] is the sum of all 2PI diagrams. 

In the absence of external sources the equation of motion for the sigma mean-field a", obtained by variation of the 
effective action | [26l l with respect to a", is 

and for the full quark propagator S"^, obtained by variation with respect to S"^, 

^=0. (29) 

The proper self energy of the quarks is 

E'"'(x,y;S) = Sg^(x,y)-i - S"\x,yr' . (30) 
Inserting the self energy (|30] l into the effective potential ( |26l l and neglecting constant terms gives 

T[cr,S] = SdM - iJrhxS-^ - /TrES + F2 S] . (31) 
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With SH/SS = the variation (|29) reads 



.,X«fc(;,,y) = _il4M^. (32) 



This means that the equation of motion for the full quark propagator S"'' | [29l l is equivalent to equation | [30l l, where 
the self energy is given by the expression | |32)| . From equation l|30] | we obtain the Schwinger-Dyson equation 

So"iS-2:S = l, (33) 

which is in explicit terms 

{id -mf)S"\x,y) -i j^d^zi:'''{x,z)S'\z,y) = i6f{x-y) . (34) 
Since implicit dependencies are not varied, the equation of motion for the sigma mean-field is 



(35) 



To solve the equation of motion for the quark propagator | |34t and for the sigma mean-field | |35] l we need the explicit 
form of the self energy and, thus, with expression | |32l l the explicit form of T2 S\ . 



A. The explicit form of T2 Wi S] and the self energy 

Since all graphs with more than one mean-field insertion are necessarily two-particle reducible, they are not in- 
cluded in T2 W> S\ . A single mean-field insertion represents them all, and we have to calculate only one diagram 
within the closed time path formalism, see figure [S] For the inclusion of such types of diagrams in the <J>-functional 
approach see [62]. Note, that this is exact within the mean-field approximation for the sigma field, because there 
are no quark self-interactions in the theory that could contribute to the two-particle irreducible effective action. It 
consists of one graph with a -|— and one with a — vertex. The corresponding Feynman rules are taken from 13411 . 

The diagram is 

^2[o-,S]=g[ d^xtr(S++(x,x)c7+(x) + S""(x,x)(7-(x)), (36) 



C 

with the trace operation tr = Eflavor Eoirac- Then, the self energy from 
reads 

i:'"{x,y) = -igSf{x-yy{x). (37) 

B. The coupled equations of motion 

With the explicit form of the self energy | [37l l the Schwinger-Dyson equation 

{id-mf)S''\x,y)-ga\x)S''\x,y) = i5f{x-y), (38) Figure 5. The only diagram for 

and the field equation for the sigma mean-field 




da" 



gtrS™(x,x) (39) 



are a coupled set of equations. In principle, we have to solve | [38t and put the solution for the full propagator S"^ 
into | [39] l. Due to the space-time dependence of cr" in < |38l l it is generally nontrivial to find the solution for S"'', which 
is exact for the given form of T2 l l36l l. We, therefore, have to approximate the full propagator. This is a crucial aspect 
because only the full approach of the two-particle irreducible effective action is a conserving, selfconsistent and 
thermodynamically consistent approximation to the exact quantum field theory |31, 32, 34]. For exact solutions of 
coupled propagator and mean-field dynamics for some model systems, see e.g. I63.-65.1 . 
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We split the mean field into one component cTq that has a slow variation compared to S"^ and a fluctuation part 
Sa", which we assume to be small. We will later ciisclose the actual meaning of this splitting 

cr''{x)=cr^{x)+Scr''{x). (40) 

We also expand the full propagator around the thermal propagator 

S"\x,y) = st{x,y) + SS''\x,y) + S^S"\x,y) . (41) 

Then for the various orders of the expansion the Schwinger-Dyson equation reads 



O{0) 
0{1) 
0{2) 



{id - mf)St{x,y) - ga',{x)St{x,y) = iSf{x - y) (42a) 
{id - mf )SS"\x,y) - ga^{x)SS''\x,y) - gScr"{x)St{x,y) = (42b) 
{id - mf)3^S"\x,y) - ga^{x)3^S'\x,y) - gSa"{x)3S"\x,y) =0. (42c) 



From Il42al l we see that the cTq part of the sigma field generates the mass of the quarks dynamically m = ntjr + gCQ. 
As already noted the idea is that the x-dependence of ctq is weak compared to the x-dependence of the propagator S. 
We identify the spatial and temporal variation of sigma with the corresponding variation of the local temperature in 
the fluid dynamic description of the quarks and antiquarks. In that spirit the solution of equation (|42a|l is given by 

zS++(p) = (Tip{x)ip^{x')) = {p + m)( \. +liTinY{\v^\)5{v'^ - m^)] forf,i'onC+, (43a) 

XP ni + le J 

iS+^{p) = ±{ip''{x')xp{x)) =2in{pi + m){nF{\p°\) -&{-p°))S{p^ -m^) for t onC+,t' onC- , (43b) 

/S"+(p) = {ipix)ip^{x')) =2iTz{p + m){nT,i\p^\) -@{p^))Sip^ -m^) for t' on C+,t on C- , (43c) 

iS—{p) = {Taip{x)ip\x')) = {v + m) ( = \ -+2inn^{\p^\)S{p^-m^)) fort,t'onC-. (43d) 

\ p fn^ — 1£ J 

From l l42bt 



SS'\x,y) = -ig j^d^zSli{x,z)5a'{z)St{^,y), (44) 



and from 



6H'\x,y) = -g" / d^zd^z'Sli{x,z')5a'{z')Sli{z',z)5a'{z)Si^{z,y) . (45) 



The approximated propagator i HTt rewritten in center and relative variables 

5a = ^{5a+ + 5a-) 
AScr = Sa+ - 5a- (46b) 



5a=^{5a+ + 5a-) (46a) 



yields 

frS''\x,y) = irSt{x,y) + ig rd^y5a{y) (s<{x - y)S>{y -x)-S>{x- y)S<{y - x) 



- '-g jd^yA5a{y) (s< {x - y)S> {y - x) + S> {x - y)S< (y - x) j (47) 

In the last step we identified the same damping kernel as in dlSt . A term similar to the noise kernel in ( [T6l l vanished 
by taking A5a = 0. In order to recover the noise kernel we need to calculate the effective action r[(j, S] explicitly 
from ||26] l with the approximations of the propagator i HTI l and | |44H45| |. It is 

Y[a,S] = SciM +gtr S++(x,x)At7(x) 

H4i,AArr('rU/Tl'i/W c;<l'r - 1/Vc;>l'i/ - r\ - ^> ( r - m'\^<(m - r\\ 

(48) 



+ /d'^ fVyA5a{x)5-a{y) [s< (x - y)S>{y - x) - S> (x - y)S<(y - x) 
/d'^ jd\jA5a{x)A5a{y) (s<{x-y)S>{y-x) + S>{x-y)S<{y-x) 
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Figure 6. The complete time path contour with the real time paths C+ and C_ and the imaginary time path C3. 



We can readily identify the same damping ( flSl l and noise kernel i fTSl l as we found in the influence functional 
approach. 

The equation of motion for the sigma mean-field obtained by varying T[a, S] with respect to Au is 

- ^^f^ = S++ix,x) + D(x) + ax) , (49) 

where we introduced the same stochastic field as was discussed in section IIII Al Note, that within the 2PI effective 
action we obtain the standard mean-field result as the first term on the right hand side of equation | |49t , which is 
formally of order g, while the damping and the noise kernel are of order g^. Therefore, the standard mean-field 
result is the lowest order contribution. 



C. The thermodynamic quantities of the quark fluid 

Using the 2PI effective action we can include the local equilibrium properties of the quark fluid. In | |48l l the calcu- 
lations were performed along the real-time contour from figure (6] For the calculation of equilibrium properties of a 
thermodynamic system we cannot, however, neglect the imaginary- time path C3. In the imaginary- time formalism 
the thermodynamic potential can be evaluated in a diagrammatic expansion [66]. The perturbative expansion of the 
thermodynamic potential in real-time formalism is more difficult. The Bogoliubov assumption of an uncorrelated 
initial state [67] leads to the factorization of the generating functional [68]. For the derivation of the Green's functions 
and the dynamics of the system the imaginary-time path C3 of the contour can be neglected. To obtain the correct 
equilibrium properties of the coupled system it must be included since it contributes to the pressure 

p = ^(InZc+lnZs) = -T[a, S] + ^ In Z3. (50) 

In equilibrium with T[a, S] \^a=0 ~ the full pressure is given by the imaginary- time path C3 of the full contour, 
see figure (6] Due to the Bogoliubov initial conditions we are left only with the one-loop effective potential. More 
advanced techniques are required to set up a consistent real-time perturbation expansion for equilibrium properties 
ll69ll . Since we restricted the model to the mean-field dynamics higher loop corrections to the pressure associated with 
propagators of the sigma field are discarded in the entire setting. This issue was briefly mentioned in the beginning 
of this section. Here, it assures that we capture the full equilibrium properties by the mean-field pressure calculated 
in equation 



V. THE EQUATION OF MOTION FOR THE SIGMA FIELD 



We now turn to the explicit calculation of the terms in the equation of motion for the sigma field | |49] |. On the right 
hand side, it includes the lowest order contribution, a damping term and the correlation of the noise. 
The free fermionic propagators on the real-time contour are given by equation ( i43l l. 
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In coordinate space the propagators obey the relations 

S++(x-y) = S-+{x-y)@{xo-ijo) + S+-{x-y)@{ijo-xo) (51a) 
S~~{x-y) = S+^{x-y)&{xo-yo) + S~^ix-y)&{yo-xo). (51b) 



A. Lowest order 



For the lowest order contribution we calculate the first term on the right hand side of equation (|49] l with the 
thermal part of the free quark propagator | |43al l 

gtiS++{x,x)=ig I Al_^2inn^{\p^\){YPi^ + m^(x))5{ip^-m^^) 

= -gPs{x) 

with the d5mamically generated quark mass mq[x) = gc{x) and the degeneracy factor dq = 12 from the trace over 

flavor, color and spin. The energy of the quarks is Ep = + and ps is the one-loop scalar density We see that 

to this order the equation of motion does not include any terms leading to damping and noise. It is the same classical 
equation of motion that was used in previous chiral fluid dynamic models. 



B. The damping kernel 

The explicit form of the damping kernel D(x) is given in equation ((TSl l. For its evaluation we define the following 
quantity 

Mix-y)=tr{S<{x-y)S>{y-x)-S>{x-y)S<{y-x)) . (53) 

Its Fourier transform can be evaluated, tracing over flavor, color and spin, and performing the integration over p^ 
by the use of the delta functions from | |43bb and | |43cl |. 

We sort the various scattering processes according to their energy balance, given by the delta functions, in order 
to make the physical processes obvious. We obtain 

M{cv,k) = -^ [ dV^^x 
47r^ J tptk+p 

X (^-{Ep + Efc+p)£p + - k ■ p)x 

X |j(a; + Ep + E;t+p)(«F(E/c+p)«F(Ep) - (1 - nF(Ep))(l - nF(E/c+p))) 

+S{cv -Ep- Efc+p)((l - nF(Ep))(l - nf{E^+p)) - n^iE^+p)n^{Ep))} (54) 
+ ((-Ep + Efc+p)Ep+2m^-k-p)x 

X \^S{a> + Ep- E^.+p)((l - nj:{Ep))nj:{E,,^p) - nf{Ep){l - n^iE^+j,))) 
+S{a} -Ep + E,t+p)(«F(Ep)(l - n^{Ek+p)) - (1 - nf{Ep))nf{E,,+p))^ 

Here, one sees the antisymmetric property Ai{—cv,\i) = — A1(a;,k). The structure of A1(a;,k) is the following. It 
describes the difference between a gain and a loss term. For example, the term np ( E)t+p ) n-p{Ep) is the probability for a 
quark-antiquark pair to form a sigma mode qq — > a, and the term (1 — Mp(Ep))(l — «p(Ej.+p)) is the statistical weight 
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of the decay of a sigma mode to a quark-antiquark pair a — t- qq. The mixed terms (1 — nF(Ep))wF(Ei:+p) describe the 
scattering of a quark (antiquark) off a sigma mode to form an anti quark (quark). For each delta function, the ratio of 
the loss to the gain term is 

^=exp(f). (55) 

^ gain \ 1 ' 

This is the detailed balance relation for the thermal quarks and antiquarks. For fixed quark masses the delta functions 
constrain the allowed scattering processes tTOll . 
The damping kernel is then 



Dix)=ig^ d'yM{x-y)Sa{y) 



.2 f d^fc 



(56) 



I d y I j^^exp{-ik{x-y))Micv,k)3a{y) , 



where we take the spatial Fourier transform 

Sd-iyo,^) = J d^yexp{-ik ■ y)Sd-{yo,y) . (57) 

To simplify further calculations we make the substitution yo = xq — t and assume that the initial time xq — yo — ?• — oo, 
such that T — 7- 00. Then 

/d^k f dco 

■^2^exp(zk-x)y^ dr J j;^exp{-ia;T)M{co,'k)6a{xo-T,k) . (58) 

We see that the sigma mean-field now depends on the history t < xq. It is known that an instantaneous approxima- 
tion to this time dependence is too constraining because the dissipative terms vanish [ 16,il7i] . We assume that the 
effect of the past can be described by harmonic oscillations around a constant value 

o-{t - T,k) ~ a{t) cos(EfcT) + b{t) sin(EfcT) . (59) 



This is the linear harmonic approximation which is also used in 119L 12211 . Here, we obtain the coefficients from the 
requirements at T = 

^(t-T,k)|^=o = ^o(i.k) ^ ait) = ao{t,k) (60) 
3r 1-° = -^ ^ bit) = --d,,it,k). (61) 



Then, we see that 



1 

a-(f- T,k) ~ (7o(f,k)cos(E,tT) - —dt(7{t,'k) sin(E;tT) 

Ej: (62) 



We now see the meaning of the splitting of u(x) that was done in ( |40] | and 



1 

Sd-{t,k) = cro(t,k)(cos(E/cT) -1) - — 3f^(f,k)sin(E;tT). (63) 



The first term gives a mass shift for the sigma field, which is only a correction to the leading order result 
We assume that this correction is small as cos(E)fT) — 1 0. Then, we can replace the fluctuation Sa{t,k.) by the 
sine-modulated time derivative of the full field. With this we can evaluate the integral over the history and obtain 
quantities that are local in time. Such an approximation will later be used for the derivation of the noise correlator, 
too. Writing 

1 

3cr{t - T,k) = -— — (exp(fEjtT) - exp{-iE^T))dtcr{t,k) (64) 
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and using the relation 



I 



dTexp(i(w — E)t) = iV — 
CO ■ 



Ti3{co — E) 



(65) 



we arrive at 



D{x) 



; iV 

2£;c V Eyt - 



do; 



M{cv,k)y 



(2^ 



exp(z'k 



(271) 

TzS{Ej- — co) — iV- 
^M{co,k)^S{cv-E,)dta{t,k). 

(27T) 



-co 



ttS{—co — Eif))dta{t,k) 



(66) 



In the final step, we used that the principle integral terms cancel by applying the antisymmetry of A4{co,k). We 
obtain for the damping kernel 



D(x) 



d^k 1 

exp(zk ■ x) —M {Ek, k)dta{t, k) . 



(67) 



A4{L0,k) contains the on-shell reaction rate of the processes given in equation ( |54l l. They lead to the dissipative 
part of the equation of motion. 

In a perturbative expansion the damping term appears first at next-to-leading order g-^ as one can immediately 
read off from (|67t . Also the thermal mass correction, which we have neglected is of this order g^. However, it is in 
fact only a correction to the mass of the sigma meson, which gets contributions from leading order | |52l l, i.e. from the 
standard mean-field contribution, and even more from the sigma field potential (|2]l. 

Being interested in the long-range oscillations of sigma we calculate the damping coefficient rj for the zero mode, 
k = 0, of the sigma mean-field and approximate (Ejt, k) ^ A4 {nia, 0) . Then for ma > 2mq only the process a ^ qq 
and the reverse reaction qq ^ a are kinematically possible. We find 



(m2-E2^ 



Mima,0) = J d'p—^^{l-2n-p{Ep))S{ma-2Ep) 



With the same approximation Ej. 

Dfx) : 



the damping kernel becomes 



m; 



3/2 



(68) 



d'^A: 1 

2 exp(zk • x)^^—M{mff,0)dtd-{t,k) 

3/2 



(27T 



(69) 



With the equation of motion | |49t the damping coefficient can be identified as 



V 



TZ \ 



2m 



2 



1 

nicr 



4 



3/2 



(70) 



Its temperature dependence is shown in figure [7] for the three different phase transition scenarios. The value of the 
sigma field in ntq = ga and the sigma mass nicr are taken to be the equilibrium values at the given temperature. 
Since rj ck g^ it is larger in a first order scenario than in a scenario with a critical point. This issue appears because 
we work at f( = and tune the strength of the phase transition by different values of the coupling g. In the linear 
sigma model with constituent quarks a realistic constituent quark mass is obtained for g = 3.3. Though the linear 
sigma model does not include confinement the damping coefficient for the zero mode of the sigma field obtained 
from the interaction with the quarks vanishes below the phase transition. This gives a realistic description at low 
temperatures. The reason is that at high temperatures the (mostly dynamically generated) quark mass is small and 
therefore mu > 2mq is satisfied. Hence the reactions a ^ qq can take place. With lower temperatures niq rises and 
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Figure 7. Temperature dependence of the damping coefficient t] for a the different couphngs g = 5.5, g = 3.63 and g = 3.3, which 
correspond to scenarios with a first order phase transition, a critical point and a crossover. 



at some point the reactions, which cause damping and noise in our model are kinematically forbidden. Physically, 
we expect that at low temperatures the decay and formation processes a ^ In become important since the pions as 
quasi-Goldstone bosons of chiral symmetry breaking become very light. In the present approach we have neglected 
the pions. This will be improved in the future. In that context it is interesting to note that our values for the damping 
are very large compared to the ones deduced from the linear sigma model without quarks Ii22,i71il . 



C. Correlation of the noise fields 

For the derivation of the correlation of the noise fields we perform the same steps for the noise kernel ( [T6t as for 
A4{x — y). For the Fourier transform of ( [T6l l we find analogously 



p'^k+p 

{-{Ep + Ek+p)Ep + 2ot^ - k ■ p) ) 



X + Ep + E;c+p)(wF(E;t+p)«F(Ep) + (1 " «F(Ep))(l " nF(Efc+p))) 

+3{u} -Ep- E)c+p)((l - «F(Ep))(l - nF(Efc+p)) + wf(E)c+p)wf(Ep))} 
+ {{-Ep + Ek^p)Ep+2m] - k ■ p)x 

X |j(a; + Ep- Et+p)((l - nF(Ep))«F(E/c+p) + nF(Ep)(l - nF(E;t+p))) 
+5{w -Ep + Efc+p)(«F(Ep)(l - «F(Efc+p)) + (1 - «f(Ep))wf(E, 



(71) 



Since the noise term has the same microscopic origin as the damping term it is not surprising that the structure is 
very similar to ( |54l l. Especially, we find that M{to, k) is proportional to the sum of the loss and the gain term of the 
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same scattering processes. The variance of the noise fields is 



d'^k (72) 
iM{co,'k) exp{-ico{t - t')) exp(zk • (x - x')) , 



where the average ()(^ is taken with respect to the Gauss distribution ( |20t . With the approximation 7V(a7, k) ~ 
N{mcr,Q) the integral over k can be evaluated to yield a spatial delta-function, which is replaced by the inverse 
volume term. Then 



{at)m)^=\ J ^Af{m,,0)expi-ico{t-t')) 



(73) 



With 



24{Ep) - 2n^{Ep) + 1 = (1 - 2nF(Ep)) coth ( ^ ) (74) 



the relation between Af{ma, 0) and M (nia, 0) is 

Afim^,0) = A^(m,.,0)coth (^) . (75) 

With the explicit form of A4 {ma, 0) | [68l l and the damping coefficient l lTOl l, we finally find 

{man)?= ^Kt-t')m,rjcoth[^) . (76) 

The approximation k = and, thus, co = nicr leads to the delta-function in the noise correlator. The noise fields 
are only correlated for equal times. This is the white-noise or Markovian approximation. In the non-Markovian 
description J\f{LO,'k) and ^A{co,'k.) have a full dependence on k. The delta-function in | [73l l is replaced by a noise 
kernel that includes the memory effects of the history of the noise fields. In addition, the damping kernel in | [69l l is 
nonlocal. Non-Markovian noises pose a difficult problem for numerical studies 1721 . 

Finally, the equation of motion for the sigma field is 

d,,d>'a+ — +gps + t]dta{x) = ^{x) (77) 

with the scalar density ps defined in from ll52t , the damping coefficient rj given in l(70l l and the correlation of the noise 
field ^ given in | [76t . 

VI. ENERGY-MOMENTUM CONSERVATION 

From the Lagrangian ^ we calculate the divergence of the total averaged energy momentum tensor. On the 
operator level we have the Dirac equation for the quark operator and the conjugate for the adjoint operator 

{i^ - gcr)q = and q{i^ + ga) = . (78) 
Then, the energy-momentum tensor for the quarks reads 

T^\x) = iq[x)Y^'q{x). (79) 

Taking the divergence yields 

dp T^" = iq^d"q + iqd" ^ q ^g^^ 
= gqqd''cr, 



17 



where we used the Dirac equation | [78l l. The energy-momentum tensor for the sigma field can easily be derived from 
the purely mesonic part of the Lagrangian (HJ 



C,r = -d^^adf'a -U{a,7i = 0). (81) 



The equation of motion for the sigma field is then found by the variational principle 



d^,df'cr+^+gqq = 0. (82) 



The divergence of the energy-momentum tensor for the sigma field is 



dfjr = ( df^df'a + ^ ) aV = -gqqd^a . (83) 



From ( |80] l and | |83l l it is clear that the sum a^iTg"' + dpT}}^ = and, thus, the total energy of the system is conserved. 
Within the full formalism of the two-particle irreducible effective action we can now take the ensemble averages of 
the calculated quantities. Since the sigma field is treated in mean-field approximation we find with | |43al l. 



and 



a„T„^''(x) = -S++(x,x)aXx) (84) 



a^rr = s++(x,x)axx) . (85) 

The total averaged energy-momentum tensor is, thus, conserved, 

'r^al = ^rW + 7^rW. (86) 

The situation is more difficult for any approximation to the full propagator due to the space-time dependence of 
the effective mass generated by the dynamic symmetry breaking. We write the energy-momentum tensor of the 
quarks in the symmetric form 

(x) = '-q{x){y''d''qix)+Yd''q{x)) - ^{d^q{x)y'' + d''q{x)jt'M^) 

-^,q{y)q{x)\x=yY -'^],q{y)q{x)\x=yY) ■ 

After a transformation to center and relative variable X = l/2{x + y) and u = x — y we see that the differentiation 
with respect to the center variable cancels and the remaining expression for the energy-momentum tensor of the 
quarks reads 

(X) = -'-d:s^-{x,u)u=o7'' - ^^!:s+-(x,«)|„=o7^ (88) 

The energy-momentum tensor of the sigma field remains the same by using the equation of motion ( |39l l defined on 
the center variable 

a,,Tr'(x) = (^a,,av + av = ^av = s++(x)a^(7(x) . (89) 

With the approximation to the full propagator l|4T)l to first order (|44)l the energy-momentum balance of the entire 
system reads 



a,.Tg%pr.(x) + a^Tf';appr.(x) = 9,, (" ^ajis^- |„=o7'' - ^3!^S+-(X,M)|„^o7 

- la^„JS+-(X,M)|„=o7" - ^d^u^s^-{x,u)i=oY) 
+ S++(X)a^<7(X) + 3S^^{X)d\aiX) , 



2 

i^,..^. s, ., (90) 
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where we identify the four contributions 

d^T^Up-iX) + a,Capp,(X) = d,T^;^{X) + dXlhi^) + h^Tll'iX) + df^STrix) (91) 
The first term evaluates to 

rj;;(x) = -^d:st^-{x,u)i=o7'' - ^a!:s+-(x,«)|„.o7'' 

It gives the energy-momentum tensor for an ideal fluid with the energy density and the pressure obtained from the 
equihbrium one-loop effective potential in mean-field approximation. This is exactly what we intend to use for the 
fluid dynamic expansion of the quark-antiquark fluid. 

In the present nonequilibrium model we find a correction, which from l|47|l is given by 

a,,^Tr(X) = D(X)a^cr(X). (93) 

With the explicit result of the damping kernel D(X) for the zero mode ( |69l , the total energy-momentum dissipation 
from the sigma field is 

di^TZppr. = i-gPs - V^t^) 9V . (94) 

It includes the dissipative part of the dynamics of the sigma mean-field. It cannot, however, account for the 
average energy transfer from the heat bath to the field given by the auxiliary noise field ^. 
What remains is the correction to the energy-momentum tensor of the quark fluid iJT^^ (X). 

In upcoming works on the numerical implementation we will investigate how well the made approximations 
conserve energy and momentum in a fully coupled dynamic system of the chiral fields and the quark fluid. 

VII. CONCLUSIONS 

In summary, we have presented a consistent nonequilibrium approach to chiral fluid dynamics, which on the 
one hand extends existing chiral fluid dynamic models by the inclusion of dissipation and fluctuations and on the 
other hand goes beyond existing studies of Langevin dynamics by putting special emphasis on the local equilibrium 
properties of the heat bath, i.e. on the back reaction of the chiral modes on the heat bath. 

The sigma field as the order parameter of chiral symmetry breaking is coupled to a fluid dynamic expansion of 
quarks and antiquarks. The interaction is given by the linear sigma model with constituent quarks, which exhibits 
dynamic chiral symmetry breaking. Due to this coupling the effective potential for the sigma field changes by the 
cooling given by the expanding quark fluid. 

We succeeded in deriving the relaxational dynamics of the sigma field from the 2PI effective action. In exist- 
ing chiral fluid dynamic models the sigma field is propagated according to a deterministic classical Euler-Lagrange 
equation of motion. The 2PI effective action includes dissipative processes and gives rise to a damping term and a 
stochastic field. We explicitly evaluated the damping coefficient and the correlation of the stochastic field in Marko- 
vian approximation for the zero mode of the sigma field. 

Although there is no confinement in the underlying theory, the damping coefficient caused by the interaction of 
the sigma field with the quarks vanishes below the phase transition temperature due to kinematic reasons. While the 
quarks gain the constituent quark mass, the sigma mass gets smaller at the first order phase transition and very small 
at a critical point. Even at a realistic coupling of ^ = 3.3 the vacuum sigma mass is larger than twice the constituent 
quark mass. 

The damping coefficient can similarily be derived in the influence functional method, where an explicit splitting of 
the system in a relevant sector, here the sigma field, and an environment, here the quarks, must a priori be assumed. 
In these terms the quark fluid acts as a locally equilibrated heat bath. In the formalism of the influence funtional, 
however, we have no control over the equilibrium properties of the quark fluid. 

In the work presented, we put special emphasis on the consistent equilibrium properties of the heat bath. This is 
the main advantage of the 2PI effective action. It is a conserving and selfconsistent approximation to the full quantum 
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theory. Besides the equation of motion for the sigma mean-field we obtain a Dyson-Schwinger equation for the real- 
time quark propagators. From the exact (for a given 2P1 effective action) solution for the quark propagator we could 
construct a conserved energy-momentum tensor. For an explicit solution to the Dyson-Schwinger equation we had 
to make further approximations. We were able to identify different terms in the divergence of the energy-momentum 
tensor of the entire system: a thermal part which coincides with the energy-momentum tensor for the classical fields 
and a correction term for both the quark and the sigma contributions to the energy-momentum balance. The correc- 
tion to energy-momentum tensor of the sigma field includes the dissipative dynamics of the mean-field. However, 
it does not account for the fluctuation energy transferred to the heat bath via the stochastic noise field. 

The presented set up gives a consistent nonequilibrium description of the coupled dynamics of the sigma field and 
the quark fluid. The entire system expands and cools. It thus describes a realistic expansion of a heavy-ion collision 
modeled by ideal fluid dynamics. Numerical results will be published in a separate work. 

Below the phase transition the zero-mode damping coefficient originating from the interaction of the sigma field 
with the quarks vanishes. It would be interesting to include the effect of higher modes and see how this leads to 
additional damping processes. However, this would go beyond the Markovian approximation and thus complicates 
future numerical studies. Additional damping processes potentially also come from the interaction of the soft modes 
of the sigma field with the hard sigma and pion modes. These processes definitely occur below the phase transition 
and assure relaxational dynamics of the sigma field. Especially the decay and formation processes a -v^ In become 
important at low temperatures where the pions are light. 

The fluid dynamic treatment of the quarks might not always be valid, e.g. in the dilute phase. Starting from the 
Dyson-Schwinger equation one can derive a Vlasov-equation for the quark-antiquark Wigner function. In i73l - [75ll it 
was solved in the collisionless approximation. It is a more challenging task to derive dissipation and noise from a 
Vlasov treatment of the quarks and antiquarks and is subject to ongoing research. 



ACKNOWLEDGEMENTS 



The authors thank Carsten Greiner and Igor Mishustin for fruitful and inspiring discussions. M.N. gratefully 
acknowledges the hospitality of the Department of Physics and Astronomy at the University of Uppsala, where 
main parts of this work were developed. M.N. acknowledges financial support from the Stiftung Polytechnische 
Gesellschaft Frankfurt am Main. This work was supported by the Hessian Initiative for Excellence (LOEWE) through 
the Helmholtz International Center for FAIR (HlC for FAIR). 



Appendix A: The calculation of the influence functional 



The influence functional method gives a reduced description of the entire system with focus on the evolution of 
the relevant variables (p. The details of the environment are eliminated by integrating out the environmental fields 
<1> in a path integral over the closed time path contour of the real-time description of finite temperature quantum 
field theory. We apply the Bogoliubov initial conditions and neglect initial correlations between the system and the 
environment. The interaction is then adiabatically turned on. As a consequence, the initial density matrix factorizes 
in system and environmental variables, pi = pf pf ■ The whole influence of the environment is then encoded in the 
influence functional, which is defined as 



exY>{iS^[cp,cp']) = /dOi| dO(pf (Oi,*;) j j DO' x ^^^^ 
X exp(/So[<l>] + fStat[<?', I'] - ;So[<l>'] - zSint[</'', <!>']) • 

Here the path integral is over all <l>(s) and ^{s)' in real time < s < with ^{ti) = Oj and <i>{t{) = 0(. 
In order to evaluate the explicit form of the influence functional we are often forced to make a perturbative expan- 
sion in the coupling between the two sectors. We expand the exponential function of Sint in equation ||A1|| 



(A2) 



exp(zSiF [<?',</'']) = I d^ij d<I>;pf(Oi,0() J V^J DO'exp(/So[0] -zSo[0'])x 

1 + i{Sint['p,^] - Sint [</'', <!>']) - liSinti'P,^] " Sint )^ + ... 
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The expansion of Siy:[(7^,o- ] becomes 

exp(zSiF[£r+,cr-]) = J Aqf j dqf j dqr J dqrpf{q+,q+;q7,q7)x 

x{\ — ig j d^x[q'^{x)q'^{x)a^{x) — q^{x)q^{x)(T^[x)) 

"^'^ I A(r(^)'?"(^)'?"(y)'?"(y)^"(^)'^"(y) 

- q+{x)q+{x)q-{y)q-{y)cr+{x)cr-{y) 

- q-{x)q-{x)q+{y)q+{y)cr-{x)cr+{y) 
+ q- {x)q- {x)q- {y)q- {y)a- {x)(7- (y)) . 

The definition of the free quark propagator is, for a,b = +, —, 

<(^,y) = {Teq"{x)q\y))o 

X exp(/So - /So [q^, q'])q'' {x)q'' (y) . 



(A3) 



(A4) 



For the explicit evaluation, we need the four-point functions that appear in the influence functional iA3i . They are 
defined in the same way as the quark propagator (|A4|| and can be obtained from the generating functional 



Z[f],rj] = Zoexp[- J^d^xd'^yi]c{x)Sc{x,y)t]c{y)] 

by subsequent differentiation with respect to the external sources fjc arid tjc . In explicit terms 

1 / iS \ / —iS \ / iS \ / —iS 



{Tq^{x)q^{x)q^{y)q^{y)) = - 



Zo V'^'7+ (^) / V'^^+ (^) / V^'/+ (y) / (y) 

S++(0)2-S++(x-y)S++(y-x) 



ij=t]=0 



(A5) 



(A6a) 



{Tq+(x)q+{x)q {y)q (y)) = — 



i6 



-iS 



-iS 



Zo \Srj+ (x) J \3fi+ (x) J V^T (y) / K^v'iy) 

S-(0)S++(0)-S+-(x-y)S-+(y-x) 



rj=t]=0 



(A6b) 



{Tq-{x)q-{x)q+{y)q+{y)) = — 



iS 



-iS 



iS 



-iS 



Zo (x) J (x) J \Sr]+ (y) / \Sfi+ (y) 
S++(0)S-(0)-S-+(x-y)S+-(y-x) 



(A6c) 



{Tq {x)q {x)q {y)q (y)) 



is 



-iS 



iS 



-iS 



Zo \Sr]- (z) J \Si]- (x) J \5r]-{y) J iv) 
S-{Qf-S-{x-y)S-{y-x). 



Then, neglecting all two-loop contributions, which cancel for S++ (0) = S (0), 



iSivW^cr-] = -\^jA'x |d4y(-S++(x-y)S++(y-x)c7+(x)c7+(y) 



+ S+-(x-y)S-+(y-x)£r+(x)£r-(y) 
+ S-+(x-y)S+-(y-xK-(xy+(y) 
- S""(x - y)S—[y - x)a-{x)a-{y)) 



(A6d) 



(A7) 
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The structure of the influence functional becomes most obvious when rewriting it in terms of the center and relative 
field variable 

a = ^{cr+ + a-) , (A8a) 

(A8b) 



We obtain 



iSi^[a,Acr] = -^^ / / d^y . 



X ^a{x)a{y){-S^^ix - y)S^^y - x) + S+-(x - y)S-+(y - x) 
+ - y)S+"(y - x) - S—{x - y)S—{y - x)) 
+ \Aa{x)-a{y){-S^^{x - y)S^+{y - x) + S+-(x - y)S-+(y - x) 



- S-+(x - y)S+-(y -x) + S-[x - y)S-[y - x)) (A9) 
:)(-S++(x - y)S++(y - x) - S+-(x - y)S-+(y - x) 
+ S-+(x - y)S+-(y - x) + S-(x - y)S-(y - x)) 



+ lAt7(x)At7(y)(-S++(x - y)S++(y - x) - S+-(x - y)S-+(y - x) 
- S"+(x - y)S+-(y - x) - (x - y)S— (y - x)) 



With the relations | |51at and jSlbb the sums of products of propagators in the brackets in | |A9t can be evaluated. We 
write = and S + = . Finally, we are left with one term that is linear and one term that is quadratic in Aa, 

iSi^[a,Aa] = V /d^x T VyA(7(x)^(y)(S<(x - y)S>(y - x) - S> (x - y)S<(y - x)) 

■' (Aio) 

+ /d'^ Jd'yAa{x)Ai7{y){S<{x - y)S>{y - x) + S> (x - y)S<(y - x)) . 
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